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K^ ■ Abstract : We calculate the lowest translationally invariant levels of the Z3- and 24- 



symmetrical chiral Potts quantum chains, using numerical diagonalization of the hamil- 
tonian for A^ < 12 and A^ < 10 sites, respectively, and extrapolating A^ — > 00. In the 
high-temperature massive phase we find that the pattern of the low-lying zero momen- 
tum levels can be explained assuming the existence of n — 1 particles carrying 2„-charges 
Q^ I Q = l, . . . , n — 1 (mass mq), and their scattering states. In the superintegrable case the 

^^ ' masses of the n — 1 particles become proportional to their respective charges: mq = Qmi. 

Exponential convergence in A^ is observed for the single particle gaps, while power conver- 
Qh. gence is seen for the scattering levels. We also verify that qualitatively the same pattern 

(— I I appears for the self-dual and integrable cases. For general 2„ we show that the energy- 

momentum relations of the particles show a parity non-conservation asymmetry which 
for very high temperatures is exclusive due to the presence of a macroscopic momentum 
rS \ Pm = (1 — 2(5/'ra)/</', where cj) is the chiral angle and Q is the 2„-charge of the respective 

C^ ' particle. 
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1 Introduction 

The 23-chiral Potts model has been introduced in 1981 by Ostlund and Huse in order 
to describe incommensurate phases of physisorbed systems P, 0], e.g. monolayer krypton 
on a graphite surface |^] . The phase structure of several versions of the model has been 
studied by various methods: Mean field, Monte-Carlo, renormalization group, transfer- 
matrix partial diagonalization and finite-size scaling of quantum chains ^, g, y, 0] . In the 
following we shall focus mostly on the quantum chain version of the model. 

In 1981-82 quantum chain hamiltonians for the chiral Potts model were obtained H, M 



via the r-continuum limit |10, 11 1 from the Ostlund-Huse two-dimensional model. These 
were not self-dual, so that the location of the critical manifold was possible only by finite- 
size scaling [^, |l2|. In 1983 Howes, Kadanoff and DenNijs |jl^ introduced a self-dual 
23-symmetrical chiral quantum chain, which, however, does not correspond to a two- 
dimensional model with positive Boltzmann weigths. Therefore in the self-dual model the 
immediate connection to realistic physisorbed systems is lost, but its peculiar mathemat- 
ical structure has attracted much interest: Howes et al. found that the lowest gap of the 
self-dual model is linear in the inverse temperature, and subsequently in ref. ||l^ it was 
shown that the model fulfils the Dolan-Grady integrability conditions [|l5| ]. (now usually 
called "superintegrability" [|l6|). A whole series of 2„-symmetrical quantum chains has 
been defined, which satisfy the superintegrability conditions [0]. The Dolan-Grady inte- 
grability conditions have been shown to be equivalent to the Onsager algebra [|^, |l^, 1£ ] , 
which entails that all eigenvalues of the hamiltonian have a simple Ising-type form. The 
well-known Ising quantum chain in a transverse field is the 22-version of the superinte- 
grable chiral Potts models. 

Much work has been done in the past few years in order to obtain analytic expressions 
for the complete spectrum of the 23-superintegrable model |16, 20 1. In the next Section 



we shall quote some of these results. The aim of the present paper is not to add to the 
analytic calculations of the special superintegrable case, but rather to use numerical finite- 
size analysis in order to investigate, how the integrable model is embedded in the more 
general versions of the chiral 23- and Z4-models. We shall study the low-lying levels of 
the excitation spectrum and investigate the possibility of a particle interpretation. 

In the neighbourhood of conformal points of isotropic theories, very simple particle 



patterns have been found by Zamolodchikov's perturbation expansion |21], which is ap- 
plicable to the non-chiral limit of the chiral Potts-models. We shall follow these particle 
patterns by varying the chiral angles, and try to find out which special properties of the 
spectrum lead, for particular parameter values, to superintegrability. Since for this pur- 
pose we have to study the spectrum through a large probably non-integrable parameter 
range, at present there is no alternative to numerical methods. It turns out that even in 
those cases which can be solved analytically, some basic features can be discovered rather 
easily through a numerical calculation, because the exact formulae are quite involved. 

The plan of this paper is as follows: In Sec. 2 we start with the basic definitions of 
the chiral Potts model for general 2„-symmetry, which then are specialized to the 23- 
and 24-cases. Sec. 3 collects some basic analytic results which are available for the generic 
2„-symmetric model. In Sec. 4 we give our detailed finite-size numerical results which 
confirm the two-particle interpretation of the low-lying spectrum of the self-dual version 
of the Z3-model. Sec. 5 discusses the analogous results for the self-dual version of the 24- 
model and shows that the low-lying spectrum is well described in terms of three elementary 
particles. While up to this point, we consider translationally invariant states only, in Sec. 6 



we present observations on the energy-momentum dispersion relations of the elementary 
particles and the effects of parity-violation on the spectrum. Finally, Sec. 7 collects our 
conclusions. 

2 The chiral 2„-Potts quantum chain 



The chiral 2„-symmetric Potts spin quantum chain [14| with N sites is defined by the 
hamiltonian 

N n-1 

H(n) = _Y^Y^^^^^k^Xa, r'; T]-^} . (1) 

j=l k=l 

Here aj and Tj are n x n-matrices acting at site j, which satisfy the relations 

aj Ff = Tj> aj J^^^' , aJ = r^ = l, a; = exp(27ri/n). (2) 

A convenient representation for a and F in terms of diagonal and lowering matrices, 
respectively, is given by 

( 0- )z,m = (5^,mCJ'~^ {J)i,m = h+\,m (mod n). (3) 

We shall assume periodic boundary conditions: rjv+i = Fi. 

The model contains 2n — 1 parameters: the real inverse temperature A and the complex 
constants a^ and a^. We shall only consider the case of H^"-' being hermitian: a^ = a^_fc 
and ak = c^n-k- ^^^ n > 2 and generic complex a^ and ak, the Perron- Frobenius theorem 
does not apply and the spectrum of H^""' may show ground-state level crossings. H^""' 
commutes with the 2„ charge operator Q = Jlfci '^j- The eigenvalues of Q have the form 
exp(27riQ/n) with Q integer. We shall refer to the n charge sectors of the spectrum of 
//("■) by Q = 0, . . . , n — 1, respectively. Parity is not a good quantum number, but //*■"■•* is 
translational invariant, so that each eigenstate of H^"^' can be labelled by its momentum 
eigenvalue p. 

For ak = ak the model is self-dual with respect to the reflection A — > A^^. If we choose 







Oik = Oik = I — icot{-Kk/n) (4) 



then the model is " superintegrable" [^2[! i-^- it fulfils the Dolan-Grady [15| conditions and 
therefore the A-dependence of all eigenvalues E{X) of H^"^' has the special Ising-like form 



m, H: 



E{\) = a + b\ + ^'imj^l + 2\cosdj + \'^. (5) 

Here a, h and 9j are real numbers. The rrij take the values rrij = —Sj, —Sj + 1, . . . ,Sj, 
where Sj is a finite integer. For details concerning formula (|5|), see flq , P2| . 
Albertini et al. [22, ^] have obtained a spin chain of the form (|I|) with 

ak = e^(2fc/"-i)'^/ sm{7rk/n) , ak = e*(2fc/"-i)v/ sin(^A;/n) (6) 

cos 99 = Acos0 (7) 

as a limiting case of an integrable two-dimensional lattice model. The Boltzmann weights 
of their two-dimensional model do not have the usual difference property |^3| and satisfy a 
new type of Yang-Baxter equations which involve spectral parameters defined on Riemann 



surfaces of higher genus. So a quantum chain (^ with coefficients (^), (^) is integrable. 
The superintegrable case is contained in (P) for ip = (j) = 7r/2. 

A number of analytic results is available for the superintegrable case. Some of these 
will be reviewed in the next Section. Starting 1987, the Stony-Brook-group has published 



a series of papers jlC, 22, 23 1 which give analytic calculations of the spectrum of the 23- 
superintegrable quantum chain. Recently, the completeness of the analytic expressions for 
the levels of the superintegrable 23-model has been shown by Dasmahapatra et al. [p^ . 
In the non-chiral limiting case 99 = (/> = of the hamiltonian (1^) with coefficients (0) 
we obtain the lattice version of the parafermionic 2„-symmetrical Fateev-Zamolodchikov- 
models WAn-^i [p5|. These quantum chains are self-dual. At the self-dual point A = 1 
and for A^ — > 00 they show an extended conformal symmetry with central charge c = 
2(n — l)/(n -|- 2). Exact solutions of these W^n-i-hamiltonian chains p6|, p7|, I 



_ B El' 29] 

have been obtained through Bethe-ansatz techniques |25, 30 1. Very recently, Cardy |31] 
has shown that a particular integrable perturbation of the critical W^n-i-models leads 
to self-dual chiral Potts models. 

Since later we shall mostly study the 23- and 24-versions of (||), we now list how 
H^^' specializes for these cases. If we keep insisting on hermiticity, H^^' depends on three 
parameters apart from a normalization [|l3[. In accordance with (0) we write: 



9 ^ 

/7(3) = -^Y. W"^"'^i + e+*^/V+ + A (e~^<^/3r .r+^^ + e+^<^/3r+r,-+i) } . (8) 



Here Oj and Fj are 3 x 3-matrices acting at site y. 

1 
w 
a;2 




(9) 



and LO = exp(27ri/3). We shall consider < A < 00; < </>, (/? < tt (for reflection properties 
of H^^' see |l^). In order not to get lost in three-dimensional diagrams, often we shall 
concentrate on two cases in which there is one relation between the three parameters: the 
integrable case (INT), where the three parameters are related by eq.(^), and the self-dual 
case (SD) where (f) = (p. For generic (p and A, the self-dual case is not known to be 
integrable. The SD case contains the non-chiral limit ip = (j) = Q, m. which we obtain the 
yV^2-model, which coincides with the Zs-standard Potts model. It has a second-order 
phase transition at A = 1, which for A^ ^ 00 is described by a conformal field theory with 
4/5 



central charge c = 4/5 [^, 33, 34 1. 

The phase diagram of the SD-chiral 23-model shows four different phases 



35|,|3|], 



see Fig. 1: for small chiral angle ip we have oscillating massive high- and low-temperature 
phases at A < 1 and A > 1, respectively, except for a small incommensurate phase interval 
around A = 1. The two incommensurate phases appear centered around A = 1 and get 
wider in A as 93 increases. 

The hermitian 24-symmetric H^'^' contains five parameters (again apart from a nor- 
malization), which we denote by A, (j), ip, (3 and (3: 



^(4) 



N 

-V2E 



-iip/2 



aj + I3aj + e 



.^^/2^3 , 



(10) 



where aj and Tj are now 4 x 4-matrices obeying (0). For simplicity, and in agreement with 
(|6|), we shall consider only the case j3 = $ = l/\/2, so that we have again two parameters 
for the non-self-dual integrable version, and two parameters A and </> = </' for the self-dual 
case. For ip = (f) = {) and /3 = /? = l/\/2 the hamiltonian ( |lO|) coincides with that of the 
Ashkin- Teller quantum chain [p^ for the special value h = 1/3 in the notations of |35]. 



From the phase diagram of the Ashkin- Teller quantum chain we know that in the case 
h = 1/3 there is just one critical point at the self-dual value A = 1. This is described by 
the W^3 rational model at c = 1, which has an extended conformal symmetry with fields 
of spin three and four. 

3 General results for the energy gaps of the 2„-chiral Potts 
model 

Apart from the exact results for the non-chiral limiting case just mentioned at the end 
of the last section, there are many exact results for the superintegrable case (||), (^ 
(equivalent to (^) with cp = (p = it/ 2): The Z2-case, which is the standard Ising quantum 
chain in a transverse field, has been solved in refs. |1^, 11, 39 1. Recently, as mentioned 
above, also the complete spectrum of the 23-superintegrable quantum chain has been 
calculated analytically [gj]. We shall first collect some important partial results which are 
available for generic TL^. 

In order to state these, we shall define all gaps with respect to the lowest Q = 0, p = 0- 
level, which we denote by So(Q = 0,p = 0). Note that ii^o(Q = 0,p = 0) is not necessarily 
the ground state of the hamiltonian. Depending on the parameters chosen, the ground 
state may be in any charge sector and may even have non-zero momentum p, since the 
model contains incommensurate phases. Nevertheless, in this section let us consider only 
gaps between -p = 0-levels. 

By AEq^i we denote the energy difference 

AEQ^i = Ei{Q,p = 0)-Eo{Q = 0,p = 0) (11) 

where Ei{Q,p = 0) is the i-th level (i = 0, 1, . . .) of the charge Q, momentum p = 0-sector. 
In the 22-(Ising)-case of (||) it is well-known that the lowest gap AEi^ has the remarkable 
property of being linear in A: 

AEi,o = 2(1 - A). (12) 



In 1983 Howes, Kadanoff and DenNijs ||13|] found that the same is true also for the self- 
dual version of the 23-chiral Potts model, eq.(|8|) at the special chiral angle ip = cj) = tx 12. 
They discovered this by calculating the high-temperature expansion up to 10th order in A, 
finding that for these special angles all higher coefficients in the expansion, starting with 
the coefficient of A'^, are zero. 

In 1 14] it was the desire to generalize (12) so as to obtain for all Z„ and all Q the 



simple formula 

AEq,o = 2Q(1 - A) (A<1, g = l,...,n-l) (13) 

which lead to choose the values (^) for the coefficients a^ and a^ and to find the superin- 
tegrability of the model. In WM the linearity of the gaps was checked through numerical 
calculations. The validity of ( |12| ) for the superintegrable 23-model to all orders of a per- 
turbation expansion in A was first shown in |23| using a recurrence formula. Later, using 



analytic methods, Albertini et al. [22, p3] for the superintegrable 23-case, and then Baxter 
poll for all superintegrable 2„, have calculated the lowest gaps of the hamiltonian eq.(||) 
in the thermodynamic limit N ^ 00. 

Using the definition (|ll|), Baxter's analytic results for the superintegrable model can 
be summarized in the form: 

AEq=i,o = 2(1 -A) (A<1) (14) 

AEq=o,i = 2n(A - 1) (A > 1) . (15) 

We now attempt to interpret the spectrum of (|l]) and the gap formula (^3|) in terms of 
particle excitations. Consider first eq.(|l3|). This would be obtained if the model contained 
only one single particle species with mass mi = 2(1 — A) and 2„-charge Q = 1. The 
lowest AQ 7^ 1-gaps would then arise from the thresholds of the scattering of Q of these 
Q = 1-particles. 

In order to check whether it is correct to describe the spectrum at if = (p = 7r/2 
this way, or whether more fundamental particles must be present, we shall study the 
zero momentum part of the spectrum as a function of the chiral angles, when these move 
away from the special superintegrable values. We shall start considering the spectrum 
down at (p = (p = 0, where we know the particle structure from the thermally perturbed 
W^n-i-iiiodels and move then towards higher chiral angles. 

The scaling regime around A = 1 of the Ww4„_i-models is known to contain n — 1 
particle species, one in each of the non-zero 2„-charge sectors Q = 1, ... ,n — 1. The 
ratios of their masses mg are [|^] : 

mq sin(7rQ/n) 

-^ = ■ / / N 1 (Q = l,...,n- 1). (16) 

77^l sm(7r/n) 

The particles with masses mq and mn-Q form particle-antiparticle-pairs. Near A = 1 the 
mass scale mi behaves as 

mi ~ (1 - A)("+2)/2n^ (;L7) 

as it follows from the conformal dimension x = 4/(n + 2) of the leading thermal operator 
of the W^n-i-model. Looking at (|T6[), we can see that in the scaling W^n-i-model 
all particles are isolated non-degenerate levels in the spectrum of their respective charge 
sector. E.g. for the particle with mass m,Q to be on the scattering threshold of Q of the 
lightest particles with mass mi, we must have Qsm{TT/n) = sin {Qn/n), which is possible 
only in the limit n -^ 00. Similarly, e.g. for n > 5, the m4-particle is below the threshold 
777-2 + 7772, and approaches this threshold only as 77 —> 00. 

We now want to follow these single particle levels in A to outside the scaling region 
and as functions of the chiral angles c/p and (p, in order to get information about what 
will happen at {/?,(^ — > |. Since near A = 1 the gaps can only be calculated numerically 
for generic (p,ip (we shall report such numerical calculations later), it is most simple to 
consider first the small- A-expansion for the lowest p = 0-gaps of the hamiltonian (1^) . Using 
the coefficients in the form (^, but not assuming (|^), to first order in A we get: 

^^Q,o = ( E «fc(l - '^'^'') ) - ^(«Q + «"-q) + • • • 

/^^^ 2sin(7rfcQ/77) . f 2k ^7TkQ\\ 2X ( 2Q \ 

' ^-^ sm [TTk/nj \ n n J sm [TiQ/n) \ n J 



\k=l 



(18) 



Eq.(0) fulfils the CP-relation 

AEQ^i{X,ip,^) = AE„ 



i(A, -¥?, ■ 



{i = 0,l,...) 



(19) 



In the sum in (^), the pairs of terms for k and n — k are equal. 

For if = and for 93 = ^ it is easy to simplify eq.(|l8|) using the trigonometrical sum 
formulae |HI|, ||] 



n-l 



y^ sin (/cvr/n) 
fe=i 
-1 
y^ sin^ (/cvr/n) 



fe=i 
n-l 



n-l 



cot — , yj sin (kirQ/n) 

^^ fc=i 



n/2, 



For (p 



and 



fe=i 

we obtain: 



AE- 



1,0 



1 / IT 

- 3 cot — 
4 V 2n 



2 cot 



37r\ 
cot — 
2nJ 



etc. 



(20) 



AE' 



2,0 



TT 

2n 


2A 

sin (vr/n) 


9 r^r^ 


3^ 2A 

1- 



h 2 cot / , X 

2n 2n sin ylTT/n) 



+ 



We identify 



"^Q 



AK 



0,0, 



(21) 



(22) 



(23) 



and conclude from ([iq ) and (^) that the mass ratio m2/mi must decrease when going 
from the scaling region A ~ 1 to A = 0. E.g. for n = 4, eq.(|l6|) gives 1112/1111 = v2, whereas 
at A = 1 from (^) and (22) we get m,2/mi ~ 1.17157. Observe, that also according (|2^), 
771-2 stays below the Q = 2-sector threshold at 2mi, and generally, the lowest levels of each 
charge sector remain isolated down to A = 0. 

We now look how the perturbation formula ( [l8| ) behaves in the superintegrable limit. 
Inserting ip = (p = tt/2, we obtain 



n-l 



n-l 



AEqa = 2 ^ sin^ (nkQ/n) - ^ sin {2TTkQ/n) cot (vrfc/n) - 2A + . . . 

fc=l k=l 

The sums can be calculated explicitly, leading to the simple result 

A^q,0 = 2(Q-A)+O(A2). 



(24) 



(25) 



Comparing ( pSJ ) to (0), we have an apparent contradiction, since in ( [25| ) the factor Q is 
not multiplying A. 

Our detailed non-perturbative numerical analysis of the spectrum for the 23- and 
24-cases (to be discussed in the next Section) shows that the single particle levels vary 
smoothly when increasing ip and 4> starting from if = (p = 0. As soon as the chiral angles 
become non-zero, the particle-antiparticle-pairs niQ and nin-Q split: rrii decreases and 
nin-i increases as p, (j) increase. Approaching ip = (j) = tt/2, m2 becomes twice as large as 
771-1, so that Sit p = (j) = '/r/2, 7712 sits just at the 7711 + rTi-i-scattering threshold. So, our 
non-degenerate perturbation theory should break down for the channel Q = 2 at p},(p > ^ 
1 13 1, and, indeed we see this in the calculation of the A^-contribution to A£'2,o- Because 



for general 2„ the explicit expression for the A^-term in (18) is too involved, here we 



specialize and give the result for the self-dual 23-case. We find 

Ai?i,o(v5, A) = 4 sin ^^ - i| cos I + A^/Cv^) + 0{X^). (26) 

where 

/(9^) = -:^(2sin^-4sin^^ + 3V3). (27) 

3 cos 99 V -J -J / 

f{ip) is smooth at c^ = ^ (we have /(§ ) = — \/3/2), but it is singular for ip — > — ^. 
Applying (|l^), from (|2|) we get AE2fl: 

AE^^oi^, A) = 4sin ^^ _ i| cos | + X^f{-ip) + 0{X^). (28) 

While formula (^) for Ai^i^o can be used in the whole range < 99 < vr (considering non- 
negative if), the corresponding expression (^) is valid only for < C/9 < ^ and diverges 
at if = (j) = ^. The A^-contribution to the analogous 24-formula will be given in Sec. 5. It 
shows a similar divergence in the expression for m2- 

Since (^) breaks down at (/? = 7r/2, we need further information in order to decide 
whether the Q = 2-particle survives at and beyond ip = 7r/2. Then this will clarify the 
question whether the single particle picture described above makes sense at the superin- 
tegrable line. 

4 Finite-size numerical calculation of the low-lying spec- 
trum of the 23-chiral quantum chain 

4.1 Convergence exponent y 

We shall now describe several detailed numerical checks of the two-particle picture from 
the spectrum in the high-temperature massive region of the 23-model. For this purpose, 
we have calculated numerically the eight lowest p = 0-levels of each charge sector of the 
hamiltonian (P) for A^ = 2, . . . , 12 sites. While in the last Section, we concentrated on the 
single particle levels, now we shall also look into the scattering states and check, whether 
the corresponding thresholds appear as expected. 

For chains of up to A = 8 sites we are able to diagonalize H^^' exactly. Using Lanczos 
diagonalization we can use up to A = 12 sites. We shall always give the results for periodic 
boundary conditions, but we have also partially checked the results with twisted boundary 
conditions. The extrapolation to A ^ 00 is done using both the Van-den-Broeck-Schwartz 



algorithm |43| and rational approximants |Q] (for details on the application to quantum 
chains and error estimation see p5[|). 

Neighbouring higher levels sometimes cross over as functions of the number of sites N , 
so that for these there is a danger of connecting wrong sequences. This can be controlled 
by calculating the leading power of the convergence, y, defined by 

AE{N) - AE{oo) = N~y + ... (29) 

and checking the smoothness of the approximants y^'- 

yN = -H^mpl!&)IHJk) (30) 



In (H) and (gg), AE{N) is a gap AEg^i calculated for N sites. 

The calculation of yN is also very useful for distinguishing between single particle levels 
and two- or more particle scattering states: For single particle levels we expect 

A^(iV) - AE{oo) = exp(-iV/e) + . . . (31) 

i.e. exponential convergence in A^ (,^ is a correlation length), so that for these the i/n 
should increase very fast with A'^. In contrast to this, for two-particle states we should 
have a power behaviour in N, more precisely, we should have lim yjv = 2. 

4.2 The spectrum for the superintegrable case if = tc/2 

We first discuss our numerical results for the 23-superintegrable case (j) = ip = it/2. Table 1 
lists results for 16 low-lying gaps for A = 0.5 in the high-temperature regime, where from 



(12) we know that mi = AEq^i^ = 1. As we have checked by repeating the calculation 
for several other values of A < 1, this value A = 0.5 is not special regarding the structure 
of the spectrum. However, there is the nice feature that because of (|l3|) for this A the gaps 
should approach integer values in the limit N ^ 00. 

In the Q = 1-sector we see an isolated lowest state at 777-1 = 2(1 — A) followed by 
a bunch of levels at 4mi. Looking into Table 2 which gives the convergence with A^ as 
parametrized by y according to (P9|), we see that indeed the mi-level (Q = l; column i = 0) 
shows exponential convergence, whereas the levels Q = 1; i = 1, 2 converge as A^~^ and 
A^~^, respectively. This indicates that the i = 1-level is due to 3-particle (mi + ttii + ttii)- 
scattering, and the i = 2- level due to 2-particle (7772 + 7772)-scattering as written in the 
bottom line in Table 1. 

The Q = 2-sector shows no isolated ground state, but starts with a bunch of levels 
around AE = 2. There is strong evidence for another bunch around AE = 5, of which 
Table 1 shows just the first level. The convergence for the further AE = 5-levels is poor 
and we do not give the numbers. A clue to the nature of the AE = 2-levels is found 
through their convergence for increasing N: Table 2 shows clearly that the level labelled 
Q = 2; 7 = has an A^-dependence which drastically differs from that of the other levels, 
indicating that this is the level of a single particle with mass r?72 = 4(1 — A) (we have 
checked the A-dependence repeating the calculation for six other values of A). So, we can 
answer the question posed at the end of the last Section: The m,2-level can be clearly seen 
in the superintegrable case, although it lies at the edge of the scattering threshold. The 
power-behaved Q = 2-states at AE = 2 should then be interpreted as ?77i + r77i-states. 
The level Q = 2; z = 5 may be mi + 7772 + 7712, which has the correct total 23-charge 
Q = 2. 

The pattern at Q = 0, AE = 3 is as expected by this picture: There should be 
77ii + r772-states and 7771 + m,i + TTii-states, both distinguished by different values of y, as 
observed in Tables 1 and 2. 

For other values of A ranging from A = 0.2 to A = 0.8 we find precisely the same 
structure. Not surprisingly, above A = 0.8 the convergence with respect to A^ is getting 
poor, since the relevant mass scale 777,1 is vanishing as A — > 1 and so at A = 0.8 it is already 
quite small. 

4.3 The spectrum off the superintegrable Hne 

After having found a simple two-particle pattern in the superintegrable case, we now check 
how this structure gets modified for 93, (/> 7^ ^. In order not to vary too many parameters. 



we consider only the INT and SD cases defined at the end of Sec. 2. For the INT-case 
exphcit formulae for the spectrum are not yet available (for some first attempts, see pql ). 
For the SD-case there may be no integrability at all. So, for the following, presently there 
is no alternative to our numerical or perturbative methods. 

We start with the SD case. In Sec. 3 we discussed already the perturbation expansion 
of the lowest gaps to order A^. This indicated that, increasing (/? from the Potts- value 
C/9 = to (/9 = 7r/2 for fixed A < 1 the mass mi decreases and m2 increases until for 
ip = 7r/2 we reach the value 1712/^,1 = 2. Above ip = 7r/2 we then should have m2/mi > 2 
which makes it difficult to isolate the Q = 2-particle in the rrii + TTii-continuum (in case 
it is still there at all). Fig. 2 shows the different patterns which we expect in the three 
charge sectors for ip < ir/2 (left hand side) and for ip > tt/2 (right). The extrapolation of 
the finite-size numbers reported in Table 3 confirms all details of these expected patterns. 
For ip = 27r/3, among the three lowest Q = 2 levels at AE = 2mi we see no exponentially 
converging level. For ip > tt/2 all thresholds are determined by mi alone, these are 3mi, 
4mi and 2m,i for Q = 0, 1, 2-sectors, respectively. 

In order to see, whether the Q = 2-particle survives in the ip > 7r/2-region, we have to 
look for more and higher levels. For small A we know from eq.(^8|) (supposing this to be 
still valid) where to search for m2 at ip > tt/2, and so we have looked slightly above the 
superintegrable line at ip = 7-7r/12, A < 0.25 among the 8 lowest levels for a fast converging 
one. Indeed, as it is shown in Table 4, the level i = 6 of the Table converges much faster 
than its neighbours. So, it is a good candidate for the ?7i2-level. However, with increasing 
A this convergence diminishes quite fast and an even/odd-N hopping takes over. This 
may be due to an increasing instability of the ?7i2-particle. 

For ip < tt/2 the threshold in the Q = 0-sector depends also on m,2, being mi + m,2. 
The left column of Table 3 shows that the numerical values for the various thresholds come 
out very well and that e.g. 7712/^-1 = 1.47755 at ip = tt /A and A = 0.5. Unfortunately, the 
determination of the corresponding values for y is quite unsafe, since for <^ 7^ ^ we have 



no exact values for lim /S.Ei available as we had in the superintegrable case from (|13|). 

So the results we obtain for the hn depend strongly on the very unprecise extrapolated 
results which are used for A£j(oo). 

4.4 The integrable case for ip> 7^ tt/2 

Turning now to the integrable case INT, our data in Table 5 show the same pattern as in 
the SD case, only the convergence is less good (for the sake of brevity we have omitted the 
finite-size data). It is not clear whether the Ai?Q=2,o-level converges exponentially at all. 
This result is somewhat surprising, because one might have expected a more clear particle 
behaviour in the integrable case. We recall that nothing is known about the integrability 
of the general SD case. 

There is not much difference between the SD-case Table 3 and the INT-case Table 5. 
E.g. observe that in the right hand part of Table 5, which gives an example for chiral angles 
above the superintegrable value, the low-lying gaps come out clearly as integer multiples 
of the single scale mi, as it is expected for threshold values if m2 > 2mi. The numerical 
precision is not very good because in the neighbourhood of the incommensurate phase the 
finite-size approximants to e.g. mi first fall with A^ and then rise again, a property which 
is not easy to handle by usual extrapolation procedures. 

Apart from the just mentioned non-monotonous behaviour of the gaps for increasing A, 
the onset of the IC-phase is not felt in our study of the p = 0-levels even at ip = cj) = 57r/6, 



A = 0.5. For still larger angles, it becomes difficult and then practically impossible to 
arrange the gaps for various A^ into plausible sequences. Anyway, it is already surprising 
that for (p = 7r/2 and higher, the low p = 0-gaps do not seem to take any notice of the 
incommensurate phase boundary, and vanish smoothly for A ^ 1. 

In the low-temperature regime A > 1 the pattern of the lowest levels is the same in 
all three charge sectors: As expected, in the limit N ^ oo the ground state is three-fold 
degenerate (one isolated level in each charge sector). Above the ground state we find for 
if < 7r/2 a gap which is A(7tt.i(A~-^) -|- ?tt,2(A~^)). Then we see a quite dense sequence of 
levels starting after a gap of 3Ami(A~^) (the same gap in all three charge sectors). If we 
would have normalized H''^' including a factor vA in the denominator of (S), the gaps 
would be just the same as those in the Q = sector at the duality reflected value of A. 
This is a generalization of Baxter's result ( p]5| ) for the superintegrable case to (/? < f . 

5 Perturbative and finite-size numerical results for the spec- 
trum of the chiral Z4-quantum chain 

In the last section we have reported our numerical evidence for a two-particle picture of 
the spectrum of the off-critical chiral 23-Potts quantum chain. In this Section we shall 
show analogous numerical evidence for a three-particle structure in the high-temperature 
regime of the chiral 24-Potts quantum chain. 

We start giving the high-temperature expansion of the lowest gaps of the self-dual 
24-quantum chain up to order A^. It reads explicitly: 

AEi^o = 2(1 + 2 sin ^~ '^ ) - 2A\/2 cos((/j/2) + X^gi^p) + ... (32) 

AE2fi = 4:V2cos{ip/2) - 2A + X^hip) + ... (33) 

AEsfi = 2(1 + 2 sin ^^ '^ ) - 2Aa/2 cos((^/2) + X'^g{-ip) + ... (34) 

where g and h are the following functions: 

1 + \/2 sin (93/2) 4sin2 99 + 5sin(/3-8\/2cos(v9/2) - 1 
^^'^^ ~ 2sin((2(^ + 7r)/4) ^ 4^/2 cos 99 cos ((^/2) ^^^' 

h{ip) = (l - ^/2cos iip/2) + 2cos2 (99/2)) (36) 

The function h[(p) is singular for (/? ^ ^, similarly, g{(p) is singular for (/? — > — ^. Of course, 
(|3J) follows from (|3^ ) by a CP-transformation, see (|l9|). At 99 = = the sectors Q = 1 
and Q = 2> are degenerate. 

At 99 = 57r/6 and A = 0, according to (|32|) the gap Ai^i^o vanishes, and above this 
value of (p the ground state of the system is in the Q = 1-sector. This feature is due to 
the particular choice (5 = l/\/2 which we made after (|lO|). If we consider more general /3, 
then this ground-state level crossing moves to (/? = 7r/2 + 2sin^^(/3/-v/2). So, for /3 > 1 the 
ground state remains in the sector Q = up to 93 = vr. 

For ip above its superintegrable value 99 = ^ (with (5 = l/v2) the scattering thresholds 
in all sectors are determined by mi = AE'i^o alone, and the "single-particle gaps" eqs.(|33|) 
and (^) move above the scattering thresholds in the sectors Q = 2 and (5 = 3, respectively. 
This is quite analogous to what happened in the Q = 2-sector of the 23-model. 
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For all A > the Q = 2-gap is larger than both the Q = 1- and Q = 3-gaps. For 
small values of A both AEi^ and AE2fi decrease with increasing ip, while AE^o grows. 
At 99 = 7r/2 the gaps are integer-spaced. 

The high-temperature-expansion is reliable for small A and the lowest levels only. So, 
for larger values of A, we have calculated the six lowest p = 0-levels of each charge sector 
for the hamiltonian (|l^) numerically. In the 24-case we are able to handle up to A^ = 10 
sites only. Correspondingly, the extrapolations are less precise than in the case of 23. 

First, we present our results for some low-lying levels in the spectrum for the super- 
integrable case of (|lO|): ip = (j) = 7r/2, j3 = $ = l/v2. Table 6 contains the results at 
A = 0.5, where from (jl^) we should have mi = 1, so that we expect the gaps for N ^ 00 
to approach simple integers. Indeed, this comes out well from our numbers. As before, 
we have checked for four other values of A that otherwise there is nothing special about 
choosing A = 0.5. 

In each of the Q 7^ 0-sectors we see one level which shows very fast convergence with 
N, and we identify this with the single-particle state of mass mg = Q. 

In the charge sector Q = we can evaluate the four lowest gaps precisely enough 
to assign them to AE = 4. There should be four different types of scattering states: 
mi+mi+mi+mi, mi+mi+m2, m2+m2 and mi+m^. While the exponents determine 
the first two types, we would have to move away from ip = cp = tx 12 in order to distinguish 
the two latter states. In the sector Q = 1 we see two excited states in addition to the 
TTii-particle state. They could be one three-particle- and one two-particle scattering-state, 
both with AE = 5. This is compatible with the expected states mi+m2+m2 and m2+ra^. 

The numerical convergence of the scattering levels and the possibility of ordering the 
levels into clear sequences in N gets worse if we consider the sectors Q = 2 and Q = 3. In 
the bottom line of Table 6 we give a few quite safe assignments. 

We now look numerically, how this particle pattern gets modified if we choose values 
of the parameters (j) 7^ vr/2, ip 7^ 7r/2 for which no analytic results are available. In Table 7 
we have collected the first few lowest energy gaps for three choices of the parameters. 

First, we look at values oi (j) = ip < tx 12. For (j) = ip = 7r/4 (which is half-way to the 
W-point) and A = 0.5 the masses of the m2- and ms-states are clearly below the scattering 
thresholds and the remaining levels can excellently be explained as multi-particle states. 
In the charge-sector Q = 2 there is one state with AE = 4.6 which seems to be unexplained 
by this pattern. However, the extrapolation A^ ^^ 00 is quite delicate and therefore the 
errors may be too small. Thus, this state could also be a 2?n,i-state. 

If we choose (j) = ip > '7r/2, the low-lying levels of the spectrum are given by mi alone. 
However, for small values of A and (j) slightly above the super-integrable line the Q = 2- 
charge-sector shows a high level that converges very fast and is not an integral multiple 
of mi. The middle part of table 7 contains the explicit values for cp = ip = Sir/ 5 and 
A = 0.10. Here, the ?TT,2-particle is still clearly visible as a high level in the Q = 2 sector. 
Unfortunately, the ?TT,3-particle is not visible in the 3?7ii-continuum. 

The right column of Table 7 shows that basically the same patterns can also be observed 
in the INT case although here the approximation is less good. 

To summarize, our numerical data supports a three-particle interpretation of the low- 
lying spectrum for the chiral 24-Potts quantum chain with ip in the range from to slightly 
above n/2. 
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6 Energy momentum relations for the single particle states 

When interpreting the excitations of the model in terms of particles, we should also look for 
their energy-momentum dispersion rule. On a lattice with A^ sites and periodic boundary 
conditions, the momentum p can take the N values 

p = -[f],... ,-1,0,1,... ,[f]. (37) 

In calculating the limit N ^ oo we use 

27r 
P = — p (-7r<P<7r) (38) 

Except for the W^n-i-case (/) = ip = the hamiltonian (||) does not conserve parity, so, 
in general the curves E{P) will not be symmetrical with respect to P — > —P. 

An expansion to first order in A gives a first orientation of the energy-momentum 
relation near A = for generic 2„. In an easy generalization of ( [T^ ) we find for Q ^ 0: 



/n-l 



AEr 



where 



fc=l / 

n-l \ 

>>,(l-..«^) 

fc=l / 

Pm = il 


rn^tp p ) 1 


sin(vrQ/n)'°^^^ ^'"^ ' " 
- 2Q/n)ct>. 



(39) 



(40) 



In the second line of (39) we have inserted the definition (|6D of ak in the A-dependent 
term. Since, to first order in A, the first term of ( |39| ) contains no P-dependence, we see 
that in the high-temperature limit, the violation of parity in the dispersion relation of 
the particle with charge Q is exclusively due to the presence of a macroscopic momentum 
Pm as given in (^0|). A particle and its antiparticle feel macroscopic momenta which are 
opposite in sign. 

In the parity conserving case if = (p = the system can be made isotropic between 
space and euclidian time rescaling the hamiltonian by a suitable A-dependent factor ^. So, 
for (fi = (/) = the particle will have the energy-momentum relation of the lattice-Klein- 
Gordon equation 



E/^ = ^J^^l + K^ (41) 

or 

E^ = ml + fK^, (42) 

where ^ is the rescaling factor, and we write 

K = 2sm{P/2), mQ = iiQ^. (43) 

The correct conformal normalization factor S, for the W^n-i-quantum chains at the 
critical point A = 1 is well-known p7|, ^]: 

^(A = 1) = n. (44) 
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Using formulae p9| ) and (]2l|), (p^), the mg and the rescahng factors for the (5 = 1- and 
Q = 2-particles of the W^n-i -chains near A = can be calculated explicitly. The mq for 
Q = 1,2 have been given in (21), (22) and ([23|), from which we find: 



vr 2A 

Ai?i,o = 2cot(— )- . (1-KV2) + .-- (45) 

2n sm (vr/n) 

TT StT 2 a 

AE„ = 2c„t (-) + 2c„t(-)- -^5-^(1 -/.V2) + ... (46) 

from which we get different rescaling factors for particles with different charge: Defining 

C = ny/X (47) 

(coinciding for A = 1 with (^)), for small A we obtain: 

?Q=i(A) = — ^C + ... (48) 

nsm(2;^) 

and 

iQ^2W= . ^^-f'l +.■■■ (49) 

^2cos2(^)-^ 

Putting numbers, we find that the main variation of the ^q over the whole range < A < 1 
is determined (for low values of n) up to 10 . . . 20% by the factor (. E.g. for 24 at A ^ 
we have 

lim^o=i/C = 0.92388... lim ^o=2/C = 0.84090 .. . (50) 

similarly for higher n. For the Ising case n = 2 we have exactly Cq=i = C for all A |4^ ]. 

For if = (j) = the two parameters ^ and mq determine the dispersion curve exactly. 
E.g. in the >V^2-case (n = 3) we find for A = 0.00, 0.25, 0.50, 0.75, 0.90: 

mQ/(l - A)^/*^ = 3.46410162, 3.58582824(1), 3.674214(2), 3.73(3), 3.8(1) (51) 

and 

C/C = 0.9428090, 0.9751702(1), 0.996697(2), 1.008(6), 1.00(1), (52) 

respectively. At A = 0.75 and A = 0.9 we observe a non-monotonous behaviour of the 
levels with A^ which gives rise to a considerable uncertainty in the extrapolation A^ ^ 00. 
At lower values of A the exponential convergence clearly sets in already below A^ = 12 
sites, so that there we may trust the extrapolation. 

For A ^ the masses determining the A^-convergence go to infinity. So for A = 
there is no A^-dependence, once a minimal value of A^, dictated by the nearest neighbour- 
interaction, is reached. This is in agreement with ([48| ) and (|49| ) and shows that the high- 
temperature expansion around A = at the same time is also a nonrelativistic expansion 
of (|4l|). For the W^2-case we have checked that the order A^-term of ( p5|) agrees with the 
second order expansion term of (|4l|): We find 

AEi^o = ^ (3 - A(2 - K^) -X'^{l-K^ + K^G) + . . .) (53) 

which fits to the form 

A^i,o = mi + (eK)V(2mi) - (eK)V(8m?) + . . . (54) 
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So, certainly to this order, the momentum dependence of the W^2-siiigle particle level 



mi is that of the Klein-Gordon equation (41) or (| 



We have no physical interpretation why for low n the ratio ^q/C varies so little with 
A (in the n = 2-case we have ^ = ("), and we find it strange for the particle interpretation 
that in general, ^q{X) is Q-dependent. 

For the 23-superintegrable case Albertini et al. |2^ have given analytic formulae and 
plots of the energy-momentum depencence of the mi -particle in the high-temperature 
range < A < 1, and for the first Q = 0-excitation in the low-temperature range A > 1. 
Here we want to give simple approximate expressions for the dispersion curves of both 
Q = 1 and Q = 2-particles, not only for the superintegrable, but also for the general 
self-dual 23-chiral model at < A < 1. 

For ip,<f) > we have no simple analytic form of the dispersion equation and so we 
have performed fits to the curves E(P) given by 

3 

E{P) =J2am COS™(P -Pm)+ Y. ^rn Sm^{P " Pm) (55) 

m=0 m=l,3 

which come out good up to A ~ 0.6. Table 8 collects some of these fitted coefficients. 

For A < 0.5 the convergence in A^ is excellent for all p if we use up to A^ = 12 sites. 
For A < 0.1 we can nicely fit the data for the momentum dispersion of the Q = 1-particle 



just using eq.(3£). For larger values of A inspection of the numerical curves and the fits in 
Table 8 shows that Pm is decreasing with increasing A, and in addition, the curves start 
being unsymmetrical with respect to P = Pm- The larger the angle 99, and the larger A 
is the more terms in the expansion (|55| ) are needed for a reasonable fit. Figs. 3 - 5 show 
momentum distributions for various values of A and if together with our fitted curves. As 
is seen from Fig. 5, for 99 = 27r/3 and A = 0.5 more terms in the expansion would have 
been needed. As we approach the incommensurate region, e.g. at if = 57r/6, A = 0.3, other 
levels get below the rrii level at | -P |> vr/2, so that the identification of the TTii-dispersion 
curve gets difficult. It may be that in this region new particle types will show up. 

7 Conclusions 

We have shown that the low-lying spectrum of the massive high-temperature regime of the 
self-dual 2„-Potts quantum chain at low chiral parameter values (p ^ 27r/3 can be described 
in terms of n— 1 particles which carry Z„-charges Q = 1, . . . ,n—l (we denote their masses 
by m-g). This is seen by studying the variation of the single-particle masses from their 
Koberle-Swieca-values at zero chirality 93 = up to and above the superintegrable value 
(p = 7r/2. For low inverse temperatures A we use a perturbation expansion. For higher 
A we concentrate on the special cases n = 3 and n = 4 and diagonalize the hamiltonian 
numerically. 

In the Z3-case the mass ratio m2/mi is shown to rise continously from m,2/mi = 1 at 
99 = to 771,2/771,1 = 2 for the superintegrable case 99 = 7r/2. In the superintegrable case the 
Q = 2-particle appears precisely at the r7ii-|-777i-scattering threshold. For 93 ^ vr, mi tends 
to zero. How far the Q = 2-particle survives for ip > tt/2 is not clear. In Table 5 we give 
evidence that the Q = 2-single particle level is still present up to A = 0.15 and if = 77r/12. 
In the superintegrable case we identify two- and three-particle scattering states through 
their different power behaviour in the chain size A^. The two single particle levels show 
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exponential convergence in N. The perturbation expansion of m2 is shown to diverge at 
order A^ for ip — > tt/2. 

At small A, the nonconservation of parity in the model results only in the appearance 
of a Q-dependent macroscopic momentum Pm = (1 — 2Q/n)(j) in the energy-momentum 
dispersion relations of the particles, which is Pm = i'/'/S in the 23-case. Our numerical 
data suggest that the average Pm decreases with increasing A, perhaps Pm ^ for A — > 1, 
but we have no simple parametrization of the effects of parity violation at large A and 
restrict ourselves to giving trigonometrical fit coefficients. 

Using the same methods we have also verified in detail that the high-temperature 
massive spectrum of the chiral Z4-Potts quantum chain can be described analogously in 
terms of three massive particles with 24-charges Q = 1, 2 and 3. In the superintegrable 
case the mass ratio equals the ratio of the charges, such that now both the Q = 2 and 
Q = 3 particles appear at scattering thresholds. As before, the scattering states can easily 
be identified by their power behaviour in the chain length N, while the single particles 
show exponential convergence with A^. 



It would be interesting to use Liischer's |48| method to obtain numerical information 



on the phase-shifts or the 5-matrix from the A^-dependence of the multi-particle states. 



Figure captions 

Fig. 1: Schematic phase diagram of the self-dual 23-chiral Potts-model as defined in 
eq.(l). 

Fig. 2: Level structure for the self-dual 23-chiral Potts-model in the high-temperature 
massive region. Left hand side: for a chiral angle (p < 'k 12 (below the superintegrable 
value) and, right hand side: (/? > 7r/2 (above the superintegrable value). In the latter case 
the thresholds are determined by m\ alone. 

Fig. 3: Energy-momentum relation of particle m\ for the 23-superintegrable case and 
various values of the inverse temperature A in the high-temperature region. The small 
crosses and triangles are finite-size values calculated for A^ = 6, . . . , 12 sites. 

Fig. 4: Same as Fig. 3, but for c/9 below the superintegrable value and both for rfi\ 
(three dashed curves) and ra2 (full curve). 

Fig. 5: Same as Fig. 3, but for several values of ip above the superintegrable value 
(towards the incommensurate region), for the Q = 1-particle mi. The shift of the minima 
of the curves to the right due to the macroscopic momentum Pm is clearly seen. 
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Tables 



A^ 


i = 1 


2 


3 


4 


5 


6 


7 


7 


3.1598983 


4.8359585 


4.4520381 


5.9303676 


5.9032240 


7.0487030 


6.3557797 


8 


3.1115649 


4.4240923 


4.2075367 


5.4868664 


5.5496714 


6.3987521 


6.0338701 


9 


3.0808417 


4.1213280 


4.0168334 


5.1252880 


5.2384091 


5.8572524 


5.7243697 


10 


3.0604137 


3.8954662 


3.8661605 


4.8297326 


4.9700825 


5.4121174 


5.4423701 


11 


3.0463166 


3.7244809 


3.7455652 


4.5867795 


4.7406004 


5.0471104 


5.1917195 


12 


3.0362806 


3.5932058 


3.6478379 


4.3856960 


4.5446740 


4.7471645 


4.9714669 


oo 


2.9998(3) 


2.986(7) 


3.01(1) 


3.02(8) 


3.05(6) 


2.9(1) 


2.95(8) 


y 


3.0(1) 


2.9(2) 


2.0(1) 


2.1(1) 


1.9(1) 


2.8(6) 


? 




3mi 


3mi 


nil +ni2 


nil +1712 


mi + 1712 


3mi 


? 



N 


i = 


1 


2 


3 


4 


7 


0.9984874 


4.5173085 


5.7905491 


6.6049711 


6.1612566 


8 


0.9994289 


4.3760023 


5.4775432 


6.0785242 


5.8456347 


9 


0.9998044 


4.2809972 


5.2349161 


5.6757753 


5.5870148 


10 


0.9999405 


4.2150750 


5.0442562 


5.3652765 


5.3749672 


11 


0.9999851 


4.1680556 


4.8924532 


5.1235743 


5.2002368 


12 


0.9999978 


4.1336962 


4.7700872 


4.9335093 


5.0552644 


oo 


1.0000000 


3.99(2) 


4.00(1) 


4.01(5) 


4.01(3) 


y 


expon. 
mi 


3.0(2) 
2mi + m2 


2.0(1) 
2m2 


2.9(2) 

2mi + m2 


1.9(1) 

27712 



N 


i = 


1 


AEQ=2,i 

2 3 


4 


5 


7 


1.9995726 


2.7901373 


4.3798245 


5.6633417 




6.1294869 


8 


1.9996761 


2.6389905 


4.0256344 


5.3106885 


6.0425834 


5.8468390 


9 


1.9998304 


2.5261089 


3.7351818 


4.9697753 


5.8323103 


5.6487379 


10 


1.9999247 


2.4399607 


3.4970152 


4.6583108 


5.5766919 


5.5065400 


11 


1.9999701 


2.3729391 


3.3008958 


4.3808694 


5.3110577 


5.4022604 


12 


1.9999893 


2.3198948 


3.1384071 


4.1365977 


5.0527923 


5.3242970 


CXD 


2.0001(1) 


2.002(7) 


2.00(5) 


2.0(1) 


2.6(5) 


4.99(1) 


y 


expon. 

7712 


2.0(1) 

27771 


1.9(1) 

27771 


1.9(2) 

2/771 


? 
? 


2.9(2) 

?77l + 27772 



Table 1: The lowest energy gaps AE'g^j, as defined in eg. (pT]) , for the Za-haniiltonian 
eq.(P) and the superintegrable case ip = (j) = 7r/2, A = 0.50. The numbers given in brackets 
indicate the estimated error in the last written digit. All calculations of the gaps have 
been performed to 12 digit accuracy and using all iV = 2, . . . , 12 sites for the extrapolation 
A^ — > oo. In order to save space, in the Tables we give less digits and omit the low-A^- 
data. Details on our determination of the exponent of convergence y, defined in (|29|), 
(|30| ) is given in Table 2 below. In the last line of each Table above, we give our particle 
interpretation of the various levels, as deduced from the result for y. 
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VN for AEQ=o^i 




UN 


for AEq= 


=l,i 


UN 


for AEq-_ 


=2,i 


N 


i = 1 


2 


3 


4 


i = 


1 


2 


i = 


1 


2 


7 


2.6413 


1.7475 


1.2814 


1.0998 


5.6670 


2.1812 


1.3354 




1.5124 


1.0720 


8 


2.6955 


1.9024 


1.3808 


1.2290 


7.2946 


2.3893 


1.4389 


2.0775 


1.5900 


1.2068 


9 


2.7348 


2.0293 


1.4594 


1.3340 


9.0948 


2.4728 


1.5230 


5.4923 


1.6503 


1.3140 


10 


2.7646 


2.1348 


1.5222 


1.4212 


11.2974 


2.5376 


1.5917 


7.7025 


1.6973 


1.4013 


11 


2.7879 


2.2232 


1.5731 


1.4947 


14.5079 


2.5883 


1.6482 


9.6994 


1.7340 


1.4733 


12 


2.8067 


2.2976 


1.6148 


1.5573 


21.8070 


2.6287 


1.6948 


11.7909 


1.7633 


1.5334 


oo 


3.0(1) 


2.9(2) 


2.0(1) 


2.1(1) 


expon. 


3.0(2) 


2.0(1) 


expon. 


2.0(1) 


1.9(1) 



Table 2: Exponents y^ of the convergence A^ ^ oo as defined in eq.(|30D, for the 23- 
superintegrable case and A = 0.5 for the lower levels given in Table 1. "expon." means 
that the fast increasing sequence of the yjy indicates exponential convergence. 



(j) = ip = 7r/4, A = 


= 0.5 


ct) = ip = 27r/3, ; 


^ = 0.5 


Q i 


A^Q,i(cx)) 


Particles 


Q i 


A^Q,i(oo) 


Particles 


1 


3.9150(2) 


mi + m2 


1 


1.75(1) 


3mi 


2 


3.85(3) 


mi + m,2 


2 


1.80(5) 


3mi 


3 


3.98(4) 


mi +7712 


3 


1.8(1) 


3mi 


4 


4.70(8) 


3mi 


1 


0.5868(2) 


= mi 


5 


3.90(5) 


mi + m,2 


1 1 


2.36(4) 


Ami 


1 


1.5834366(1) 


= mi 


2 


1.174(2) 


2mi 


1 1 


4.68(1) 


2m,2 


2 1 


1.18(2) 


2mi 


2 


2.339607(1) 


= 7712 


2 2 


1.19(5) 


2mi 


2 1 


3.168(4) 


2mi 








2 2 


3.17(5) 


2mi 









Table 3: The lowest energy gaps AEq^i (extrapolated N — > 00) of the self-dual 23- 
model (^), for A = 0.5 and two values of (p: (p = 7r/4 (below the superintegrable line), and 
ip = 37r/2 (above the superintegrable line), together with their particle interpretation. 





^Eq-_ 


--2,t iv = 77r/12, A = 
mi = 1.50335845 


0.10) 


^Eq-_ 


=2,i (^ = 77r/12, A = 
mi = 1.41131611 


0.15) 


N 


i = 


5 


6 


7 


i = 


5 


6 


7 


7 


3.03919 


3.58623 


3.89205 


7.99416 


2.86777 


3.63650 


3.95076 


7.69924 


8 


3.03222 


3.72696 


3.89348 


7.90458 


2.85847 


3.85581 


3.98085 


7.57893 


9 


3.02727 


3.65640 


3.89285 


7.83788 


2.85177 


3.75262 


3.96189 


7.49035 


10 


3.02363 


3.74332 


3.89313 


7.78699 


2.84678 


3.89547 


3.97910 


7.42310 


11 


3.02088 


3.69394 


3.89301 


7.74732 


2.84297 


3.81857 


3.96735 


7.37073 


12 


3.01875 


3.75171 


3.89306 


7.71579 


2.83999 


3.91801 


3.97795 


7.32905 



Table 4: Selected Q = 2-gaps for ip = 77r/12 (slightly above the superintegrable line) 
of the self-dual 23-niodel. The numbering of the levels (i = 0, . . .) refers to iV = 12 sites, 
for a smaller number of sites there are less levels between i = and i = 5. The values of 
77T-1 are very well determined because of fast convergence. Observe that the i = 0-levels 
converge excellently to 2r7T,i. In both cases level i = 6 converges much faster with N than 
the other levels, which may indicate that this is the m2-level. Non-monotonous behaviour 
in N is common in the neighbourhood of the IC-phase. 
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99 = 7r/4, (/) = 19.47...° 


f 


= 27r/3, (/) = 135.58...° 


A = 0.75 






A = 0.70 




Q i 


A^Q,i(cx)) 


Particles 


Q 


i 


AEQ,i(oo) 


Particles 


1 


1.933(2) 


mi + ?TT,2 





1 


1.50(1) 


3mi 


2 


1.9(2) 


mi + m2 





2 


1.53(5) 


3mi 


3 


1.8(3) 


mi + m2 


1 





0.495(2) 


= mi 


1 


0.750(1) 


= mi 


1 


1 


1.99(1) 


Ami 


1 1 


2.41(2) 


2m2 


2 





0.995(1) 


2mi 


2 


1.19(2) 


= m,2 


2 


1 


1.00(1) 


2mi 


2 1 


1.50(3) 


2mi 


2 


2 


1.04(5) 


2mi 








2 


3 


2.52(5) 


5mi 



Table 5: Lowest energy gaps AEq^i as in Table 3, but for two different choices of the 
parameters in the integrable Zs-model case, in which if and </> are related by (0). On the 
left side we give an example with chiral angles below the superintegrable case, on the right 
hand side (j),(p are taken above their superintegrable values. The pattern observed here is 
very similar to that of the self-dual case in Table 3. 



N 


i = l 


AEQ=o^i 
2 3 


4 


i = 


A^Q=l,i 

1 


2 


4 


4.4870349 


6.6401198 


6.9386000 


7.1202527 


0.9782260 


6.2907501 


8.8182343 


5 


4.2443835 


6.2079844 


6.3448347 


6.4390786 


0.9932173 


5.7339386 


7.9355672 


6 


4.1352996 


5.6907770 


5.9280172 


6.0208759 


0.9985794 


5.4429212 


7.2245463 


7 


4.0807817 


5.3191415 


5.5469396 


5.7278875 


0.9999564 


5.2809826 


6.8098311 


8 


4.0511625 


5.0477033 


5.2606178 


5.2880905 


1.0001440 


5.1860132 


6.4550409 


9 


4.0339657 


4.8460464 


4.9582447 


5.0427288 


1.0000927 


5.1277432 


6.1875836 


10 


4.0234306 


4.6936534 


4.7248803 


4.8744672 


1.0000390 


5.0905398 


5.9834045 


oo 


4.0000(2) 


3.89(2) 


3.9(9) 


3.7(4) 


1.0001(1) 


4.97(3) 


5.7(6) 


y 


4.0(2) 
4m 1 


2.1(3) 

2m2 or 

mi +m3 


2.9(4) 
2mi + m2 


1.6(1) 

2m2 or 

mi + ms 


expon. 
mi 


3.5(3) 
? 


1.8(1) 
fn2 + ms 



N 


i = 


=2,i 

1 


i = 


A^Q=3,^ 

1 


2 


4 


1.9760981 


3.6344881 


3.0090309 


5.0706873 


5.4135243 


5 


1.9885583 


3.1951141 


2.9978592 


4.4899244 


4.9297443 


6 


1.9959972 


2.9027708 


2.9980224 


4.1048393 


4.5520862 


7 


1.9989351 


2.7013088 


2.9991277 


3.8412387 


4.2636142 


8 


1.9998278 


2.5582804 


2.9997209 


3.6558602 


4.0432463 


9 


2.0000213 


2.4539476 


2.9999379 


3.5222475 


3.8733462 


10 


2.0000330 


2.3759111 


2.9999964 


3.4237046 


3.7406336 


oo 


2.00001(3) 


1.9(1) 


3.00001(6) 


2.8(2) 


3.02(4) 


y 


expon. 
m-2 


1.9(1) 
2mi 


expon. 
m3 


2.1(2) 

m,i + 77l2 


1.7(2) 
mi + m2 



Table 6: The lowest energy gaps AE'q^j, as in Table 1, but here instead of the 23-case 
now for the 24-hamiltonian eg . (|lO|) . Superintegrable case cj) = ip = tt/2 for A = 0.5. Note 
the overshooting in the approximants for ttii and m,2- 
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(f) = (f = 7r/4 




(f> = If = 37r/5 




ip = Tr/4:, (j)= 19.47...° 




A = 0.50 






A = 0.10 




A = 0.75 


Q i 


A^Q,.(oo) 


Particles 


g 


7 


AEq4(X)) 


Part. 


Q i 


A^Q,.(oo) 


Particles 


1 


5.77(2) 


mi +m3 





1 


5.1(2) 


47771 


1 


2.767(4) 


7771 + "is 


2 


5.76(8) 


mi +m3 


1 





1.2116781 


= 77ii 


2 


3.2(1) 


27711 + rn2 


3 


7.1(2) 


2m2 


1 


1 


6.1(3) 


5mi 


1 


0.8663(1) 


= 7771 


1 


2.03383(1) 


= TTT-i 


2 





2.42334(7) 


27771 


1 1 


3.50(1) 


27711 + "T-S 


1 1 


7.1958(4) 


m2 +m3 


2 


1 


2.430(6) 


27771 


2 


1.529(2) 


= 7772 


1 2 


6.98(7) 


m2 +m3 


2 


5 


3.211450(1) 


= 7772 


2 1 


1.8(1) 


27771 


2 


3.4282(1) 


= 7712 


3 





3.6349(5) 


37771 


3 


1.786(2) 


= 7713 


2 1 


4.15(6) 


2mi 










3 1 


2.52(7) 


37711 


2 2 


4.59(1) 


? 
















2 3 


7.62(7) 


2m3 
















3 


3.7501(1) 


= 7713 
















3 1 


5.458(2) 


TH-i +7n,2 

















Table 7: The left and middle parts of this Table show a selection of lowest energy 
gaps of the self-dual, but not super-integrable 24-haniiltonian for two different choices of 
the parameters (/) = ^p, with (5 = $ = l/\/2 and different A. The right part of the Table 
contains an example of the INT case. For most levels in this Table the determination of the 
convergence exponent y is very unsafe due to the considerable uncertainty in A£'q^j(cxd). 
So, as also in Tables 3 and 5, the particle content is inferred almost exclusively from 
A£'q^j(oo) alone. 



[deg] 


A 


ao 


ai 


02 


03 


61 


63 


[deg] 


45. 


0.25 


2.8745 


-0.5924 


-0.0549 


-0.0075 


0.0033 




39. 


45. 


0.50 


3.0302 


-1.1653 


-0.2355 


-0.0764 


-0.0134 




30. 


45. 


0.75 


3.3502 


-1.6702 


-0.6028 


-0.2300 


0.0864 


-0.2121 


18. 


90. 


0.25 


2.0436 


-0.5118 


-0.0412 


-0.0848 


-0.0112 




81. 


90. 


0.50 


2.2363 


-1.1331 


-0.2381 


0.0014 


-0.0188 




69. 


90. 


0.70 


2.4363 


-1.4587 


-0.4244 


-0.1190 


-0.0757 


-0.0112 


57. 


90. 


0.85 


2.6161 


-1.5982 


-0.5752 


-0.2791 


-0.1981 


-0.0745 


39. 


120. 


0.25 


1.4317 


-0.5570 


-0.0552 


-0.0022 


-0.0263 




105. 


120. 


0.50 


1.5407 


-0.8112 


-0.0909 


-0.2517 


-0.3291 




57. 


150. 


0.20 


0.7281 


-0.4696 


-0.0028 


0.0451 


0.0013 




150. 


171. 


0.10 


0.2257 


-0.2260 


-0.0128 


-0.0009 






171. 


45. 


0.50 


3.7927 


-1.1973 


-0.1501 


-0.0874 


-0.3223 




-45. 



Table 8: Coefficients of fits to the energy- momentum relation eq.(|55|) for various ip and 
A, all for the self-dual model. Except for the last line, which is for the lowest Q = 2-gap, 
all other fits are for the lowest Q = 1-gaps. The quality of the fits can be judged from 
Figs. 3, 4 and 5. In the last column we quote pm^ which we define by Pm = '/'m/3. 
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